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Half integer Josephson vortices in 0-lt-junctions, discussed theoretically and observed experimentally, spon- 
taneously appear at the point where the Josephson phase is Jl-discontinuous. The creation of arbitrary disconti- 
nuities of the Josephson phase has been demonstrated recently. 

Here we study fractional vortices formed at an arbitrary K-discontinuity, discuss their stability and possible 
ground states. The two stable states are not mirror symmetric. Furthermore, the possible ground states formed 
at two K-discontinuities separated by a distance a are investigated, and the energy and the regions of stability 
of each ground state are calculated. We also show that the ground states may strongly depend on the distance 
a between the discontinuities. There is a crossover distance a c such that for a < a c and for a > a c the ground 
states may be qualitatively different. 
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I. INTRODUCTION 

During the last years it was shown both theoretically<ii2i2ii£ 
and experimentally^i^ that one can create and study half in- 
teger Josephson vortices which carry only half of the magnetic 
flux quantum <t>o- To create such vortices one usually has to 
use a so-called long Josephson 0-71 -junction (0-7I-LJJ), i.e., a 
junction of which one part behaves as a 0-junction (positive 
critical current) and the other part as a 7t-junction (negative 
critical current). There are several approaches and technolo- 
gies which can be used to fabricate such junctions. For ex- 
ample, one can use junctions based on superconductors with 
anisotropic order paramete^iSiii, junctions with ferromag- 
netic barriesi^ii^, or even conventional junctions with a pair 
of tiny injectorsi2ii£ii&. 

In fact, all this types of Josephson junctions may be de- 
scribed by a model in which the Josephson phase is 7t- 
discontinuous at the 0-7i-boundary. In an infinite LJJ, the 
presence of a discontinuity results in an infinite energy, which 
of course cannot be allowed in nature. To save energy, the 
Josephson phase (j) bends around the discontinuity (on the 
length scale of the Josephson penetration depth Xj) so that 
(|)(±oo) = 2nn. This localizes the energy in a region of size 
~ Xj around the discontinuity and creates a local magnetic 
field oc (j) v (x) and currents °< sin(j)(x) circulating around the 
discontinuity. Circulating currents create a vortex with the to- 
tal magnetic flux equal to < t ) o/2. 

Instead of specifying the flux carried by the vortex, we just 
denote it by the total continuous change of the phase on the 
interval from x = — °° to x = +°° (not including the disconti- 
nuity), i.e., a 7t-vortex is a semifluxon carrying <E>o/2, while a 
27t-vortex is the usual integer fluxon carrying <I>o. 

Recently it was demonstrated that it is possible to create 
an arbitrary K-discontinuityi2iii The question, therefore, is: 
what kind of vortices can be formed at such an arbitrary K- 
discontinuity? As we saw above, the main reason of the frac- 
tional vortex formation is to save energy (localize the energy 
around the phase discontinuity point) by compensating the 
phase jump at the discontinuity towards an integer number 



of 2n. In the case of a +7i-discontinuity, the vortex is either 
a — 7t-vortex with (|)(+°°) — (j)(— °°) = or a +7t-vortex with 
( K+ 00 ) — <K — 00 ) = 2JC. The ±7t-vortices are mirror symmet- 
ric, i.e., one carries +<Pq/2 with currents circulating clock- 
wise, while the other carries — <t>o/2 with currents circulating 
counterclockwise. 

We can already anticipate that in the case of an arbitrary 
—K-discontinuity, at least two, generally not mirror symmetric 
vortices can exist: a +K-vortex with (|)(+°°) — (j)(— °°) = and 
a (k — 27i)-vortex with (j)(+°°) — (j)(— °°) = — 2n. Below we 
present an analysis of the possible vortex states pinned at one 
or two arbitrary discontinuities of the Josephson phase. 

The paper is organized as follows. In Sec. [H] we introduce 
the model for the case of arbitrary discontinuities. In section 
IHIIwe derive an analytical expression for the phase of various 
vortices pinned at an arbitrary discontinuity. Then in Sec. IIVI 
we study numerically the ground state of a LJJ with two K- 
discontinuities situated at a distance a from each other. Sec- 
tion|V]concludes this work. 

II. THE MODEL 

The behavior of the Josephson phase §(x) in a LJJ with 
discontinuities can be described by the following perturbed 
sine-Gordon equation^ 

tyxx-tytt - sin((j)) = acj), -y(x) + d xx (x), (1) 

where the subscripts x and t denote a partial derivative with 
respect to coordinate x and time t, respectively; 0(x) is a 
step function which is K discontinuous at the points x = x, 
(i = 1 , . . . , AO of the junction and is constant everywhere else. 
We note here, that for a ^-discontinuity the sign of the jump 
does not matter. A junction with a +jt, +7C, . . . , +7t set of dis- 
continuities physically behaves in the same way as, e.g., a 
junction with +Jt, — K, +7t, — 7t, . . . , +Jt, — Jt. In contrast, the 
sign of a K-discontinuity does matter. As we will see later, 
a junction with two (— K, — k) discontinuities may have com- 
pletely different ground states than a junction with (+K, — k) 
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discontinuities. If the discontinuity is created by injectors 14 , 
the sign of discontinuity depends on the polarity of the in- 
jected currents. 

Equation Q is written in normalized units. The coordi- 
nate is normalized to the Josephson penetration depth Xj, the 
time is normalized to the inverse plasma frequency ft)" 1 , the 
bias current density y = j / j c is normalized to the critical cur- 
rent density, and a = 1 / a/P7 is the dimensionless damping 
parameter (P c is the McCumber-Stewart parameter). 

Looking at Eq. Q one can anticipate that since (j) vv and 8 VV 
are two additive terms, the solution <j)(x) should be discontinu- 
ous at x = x;. It is not very practical to deal with discontinuous 
functions such as (])(x), and with singular functions such as 
9xc(x). Therefore, it is convenient, following Ref. ,5, to intro- 
duce a new continuous phase /j: §(x,t) = n{x,t) + Q(x). Then 
the sine-Gordon equation reads 



Mvx - Utt - sin^ + 9) = a/j, - y(x) . 



III. VORTEX AT A SINGLE k-DISCONTINUITY 



(2) 



For the investigation of the ground state of a system with 
one K-discontinuity at x\ = we use the static version of 
Eq. (0 without bias current (y = 0) 



sin^ + S), 



and with 



9(x) 



0, x<0; 
-K, x> 0, 



= -kh(*), 



(3) 



(4) 



where H(x) is a Heavy side step function. 

The results in the case of a +K discontinuity are mirror 
symmetric, i.e., /u(x,k) = — fi(x,— k). Since the phase is de- 
fined modulo 2k, without loosing generality, below we con- 
sider only < K < 2k. 

The natural boundary conditions (BCs) at x = are 



M+o) = M-0); 
m*(+o) = m*(-o). 



(5a) 
(5b) 



The BCs at x — > ±°° come from the fact that the energy 
density of the system 



A- 



[l-cos^ + 0)], 



(6) 



must vanish at x — > ±°°. This can be achieved when both (j) v — > 
and (j) — > 2Kn. Since (]) A -(±°°) = ^(±0°), (|)(— °°) = /«(— °°) 
and (j)(+o°) = — K, without loosing generality we adopt 

the following BCs 



Mt(±~) = 0; 
/i(-oo) = 0; 
<"(+°°) = k + 2jc«. 



(7a) 
(7b) 
(7c) 



First we consider a +K- vortex, e.g., n = in Eq. ( TTcl . The 
non-trivial solutions of Eq. for x < and x > are always 
given by the "tails" of a fluxon, which may be shifted along 
the fj or x axes to satisfy the BCs 4 , thus 



H{x) 



<|)o(x-xo), x<0, 

K-<|)o(-X-Xo), X>0, 



where <|>o(x) is a soliton (fluxon) solution 
<|>o(x) = 4arctane Jc . 



(8) 



(9) 



Note, that the ansatz (|8j already satisfies the BCs 0. By 
trying to satisfy the BCs (0 we arrive at the expression for xq 



xq = — In tan — > 0. 
8 

Second, we consider a (k — 27i)-vortex, e.g., n 
Eq. dTcl . A similar procedure gives 



tyo(x-xo) 

K — 271 — X — Xq) 



x<0, 
x>0, 



with (j)o again given by Eq. and 



/K 7C\ 

= - In tan ( - - - J > 0. 



(10) 
-1 in 

(11) 
(12) 



In principle one can try all other n in Eq. (TTcl . but as soon as 
|k + 27C77 1 > 47C the solution cannot be constructed at all. Here, 
formally, the value of xq, which in the general case is given by 



x 



In tan 



UK 

T 



(13) 



becomes complex. Physically this means that one cannot con- 
struct a localized energy solution out of two fluxon tails so 
that the phase changes by more than 471 when x goes from — °° 
to +00. Moreover already for |k + 27Ch| > 27C the constructed 
solution is unstable^. 

Thus, for a —K-discontinuity we have only two possible 
stable fractional vortices, given by Eqs. and (II It . which 
can localize the energy around the discontinuity. In Fig.^ we 
show both phases (j)(x) and fi(x), the magnetic field /j a (x) and 
the supercurrent sin(|)(x) for both of these solutions for the 
particular case of K = k/2. Similar pictures for semifluxons 
(k = 7C) can be seen in Ref.|5| 

In the following, a K-vortex pinned at a —K-discontinuity is 
called a direct vortex and denoted as ] , while a (k — 27t)-vortex 
pinned at a —K-discontinuity is called a complementary vortex 
and denoted as [. We intentionally avoid the word "antivor- 
tex" because when at K = 2k the direct vortex is a fluxon, the 
complementary vortex is not an antifluxon. It is just a "no 
vortex" state, /j = 2Kn. By definition, the complementary of 
a complementary vortex gives again a direct +K-vortex. In 
the general case there are only two stable vortices: a direct 
vortex and a complementary vortex. The single exception is 
the case of K = 27t«, for which there exist tree stable solu- 
tions: +27C-vortex (fluxon), constant phase state (zero phase), 
and — 27C-vortex (antifluxon). The complementary "vortex" 
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FIG. 1 : Phases ty(x) and fj(x), the magnetic field fi x (x) and the supercurrent sin(])(x) for K = 7t/2. (a) direct lt/2-vortex (b) complementary 
— 37t/2-vortex tl It . The lines denoted as "fluxon" and "antifluxon" show the corresponding curves for 27t-vortex (a) and — 27t-vortex (b) for 
comparison. 



for both ±271- vortices is a constant phase state and, vice versa, 
for the constant phase state the complementary vortex is either 
a fluxon or an antifluxon. In the majority of situations this 
can be distinguished due to the conservation of the topologi- 
cal charge. Note that in the notations such as ] J \[ the direction 
of the arrows shows the polarity of the vortex (up or down), 
while the harpoon on the left or right side indicates whether 
the vortex is direct or complementary. The notations are sum- 
marized in Tab. |I] 

The energy of the direct vortex can be easily calculated by 
integrating the energy density (|6j over the junction length with 
/j(x) given by 0, which gives 

E(k) = 16sin 2 -. (14) 
8 

The energy of a complementary vortex is obviously £(271 — 
k). Note, that this expression gives the right energy of a fluxon 
£■(271) = 8 and of a semifluxon E(n) =4(2- v^), see Ref.d 



Symbol Discontinuity Topological charge Name 



1 


— K 


+K 


direct 


I 


+K 


— K 


direct 


t 


+K 


271 — K 


complementary 


I 


— K 


K— 2tc 


complementary 


T 


±71 


+71 


semifluxon 


i 


±71 


-7t 


antisemifluxon 


t 





+271 


fluxon 


4 





-271 


antifluxon 



TABLE I: Notations for different types of vortices. 



One may use energy arguments to make some conclu- 
sions about the stability of various states*. For example, can 
"heavy" vortex emit a fluxon and turn into a "lighter" compli- 
mentary vortex? If the energy of a fluxon together with the 
energy of a complimentary vortex is larger than the energy of 
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an original vortex, this process is clearly impossible. Using 
Eq. dl4t we can immediately see that such process is forbid- 
den for vortices with the topological charge less than 2k. On 
the other hand the (K + 27c)-vortex with the topological charge 
larger than 2k may emit a fiuxon and turn into a vortex with 
the topological charge K. Indeed, using Eq. (I14> one can prove 
that 

E(k + 2k) > E(2K)+E{K)fov 0<k<2tc. 

We stress that this is necessary, but not sufficient condition 
for vortex instability. Similarly, the (k— 4k )- vortex with the 
topological charge less than —2k, may emit fiuxon and turn 
itself into a (k — 271 ) -vortex. More strict analysis^ shows that 
indeed the ^-vortices with > 2?t are unstable. 

In comparison with a semifluxon, for which K and (k — 2tc) 
vortices are mirror symmetric, in the case of an arbitrary dis- 
continuity the symmetry is broken. This may result in a num- 
ber of interesting consequences both for ground states and for 
vortex dynamics. 

We would like to note here that many authors in the context 
of unconventional superconductivity propose that the pres- 
ence of a fractional vortex (not equal to <t>o or &o/2) is the 
signature of time-reversal symmetry violation. In our case, the 
time-reversal symmetry is preserved: upon inversion of time 
a K-discontinuity becomes a — K-discontinuity, and the cor- 
responding direct and complementary vortices change their 
signs too. In our case, only the parity is violated, i.e., direct 
and complementary vortices pinned at a fixed discontinuity 
are not mirror symmetric. 



IV. GROUND STATE OF TWO VORTICES 

A. Possible states 

Now we consider two K-discontinuities (0 < K < 2tc) in an 
infinite LJJ situated at x\o = ±a/2 (at a distance a from each 
other). If both discontinuities have the same sign of K, e.g., 
(— K, — K) and 



0(x) = -K 



H 



a 

X+ 2 



H 



a 

x 

2 



-27m, (15) 



there are two possible irreducible vortex configurations: 
the symmetric ferromagnetic (FM) state ]] consisting of 
(+K, +K)-vortices, and the asymmetric antiferromagnetic 
( AFM) state 1 [ consisting of (+K, K — 27c)-vortices. If the dis- 
continuities have different sign, e.g., (+K, — k) and 



8(jc) = K 



H 



a 

X+ 2 



H 



a 

X ~2 



+ 2kh, (16) 



there are two other irreducible vortex states: the asymmetric 
FM state |"| consisting of (271 — k,k)- vortices, and the sym- 
metric AFM state J"| consisting of (— K, +K)-vortices. Below 
we consider these four ground states. 

The analytical description of the two vortex states is pos- 
sible but involves rather bulky and intuitively not clear ex- 
pressions with elliptic functions. Therefore, here we present 



only a qualitative picture of the ground states and calculate the 
necessary quantities numerically. For these numerical calcu- 
lations we have used StkJJ software^. To obtain the mag- 
netic field profiles we started the simulations at K = K with the 
phase profile approximately corresponding to the one of four 
states. Then K was changed in small steps Ak = 0.017C towards 
or towards 2tc. After changing K the program was waiting 
for the decay of all oscillations in the system and, after this, a 
snapshot of the magnetic field fJ x (x) was produced. 

For an infinite LJJ containing more than one discontinuity 
and more than one vortex the topological charge of a vortex 
is not very well defined, especially if the vortices are close 
to each other. The topological charge of a vortex is given by 
/u(+oo) — 00). This definition assumes that there are no 
other topological excitations in the vicinity of a vortex. When 
vortices are close to each other it may be difficult to sepa- 
rate them, their magnetic field and currents may overlap and 
superpose (cancel or enhance) non-linearly. Nevertheless we 
will still continue using the terms K-vortex or K — 27C-vortex, 
especially in the case of weak coupling (large and moderate 
distance between vortices), assuming that the state is con- 
structed from the corresponding single vortices at an infinite 
distance from each other and then the distance between dis- 
continuities is smoothly and quasi-statically decreased down 
to desired values. We also assume that during this decrease no 
abrupt reconfigurations occurred in the system. 

First, we consider the case of a rather large a = 3 (weak 
vortex-vortex interaction), for which the states can be con- 
structed out of single vortex states. The value of a is given in 
normalized units, i.e., in units of Xj. The numerically calcu- 
lated magnetic field profiles corresponding to the four differ- 
ent states \ [, 1 1 , H and J ] are shown in Fig. [2] 

We introduce the term complementary configuration or 
state which can be obtained from any original configuration if 
all vortices are substituted by the complementary ones. Note 
that the profiles p. L x of the complementary states [], [[, \[ and 
\ [ pinned at the same discontinuities are given by 



/4(*,k) = -ii x {x,2k-k) 



(17) 



If we discuss some property, e.g., instability, that happens with 
some state at some value of K, the same thing happens with 
the complementary state at 27C — K. Thus all that is presented 
below can be easily mapped to the complementary states using 
K — > 27t — K and Eq. ( fTTl . 

a. Asymmetric AFM state. The field profiles fj x (x) 
formed at the ( — K, — k) discontinuities and corresponding to 
the asymmetric AFM state ] [ are shown in Fig.|2ja) for several 
values of K. This state is stable for any value of < K < 2k. 
At K = or at K = 27C the state becomes unstable, the integer 
fiuxon is emitted and the state turns itself into a symmetric 
FM state. These processes can be written as 



K=0„ 



4, 



fr 



b. Symmetric FM state. The field profiles fi x (x) formed 
at the (— K, — k) discontinuities and corresponding to the sym- 
metric FM state ] ] are shown in Fig. [2b) for several values of 
K. When the vortices become "heavy" their top gets rounded. 
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FIG. 2: Magnetic field profiles fi x (x) for different values of K for the states (a) ] [, (b) 1 1 , (c) |"| and (d) J 1 with discontinuities at x = ±a/2 and 
a = 3. The function 0(x) is shown by gray line. The amplitude of 0(x) is arbitrary. 



This state is stable only for < K < k' c ' , with K c ' > 7L When 
the vortices become as "heavy" as two fiuxons, they loose sta- 
bility because they repel each other and the pinning is very 
weak. As a result one integer fluxon is emitted and the sym- 
metric FM state turns itself into an asymmetric AFM state dis- 
cussed above. The value of K c - depends on a and the K c (a) 
dependence is presented and discussed below. 

c. Asymmetric FM state. The field profiles fJ x (x) formed 
at the (+K, — k) discontinuities and corresponding to the 
asymmetric FM states |"| are shown in Fig. 0c) for several 
values of K. This state is stable for any value of < K < 271. 
At K = or at K = 2n the state becomes unstable, an integer 
fluxon is emitted and the state turns itself into a symmetric 
AFM state. These processes can be written as 

t1^J1+t or if^H + f. 

d. Symmetric AFM state. The field profiles /u x (x) 
formed at the (+K, — k) discontinuities and corresponding to 
the symmetric AFM state J ] are shown in Fig.|2Id) for several 

values of K. This state is stable only for < K < 1Q 1 , iq > jr. 

1 1 

When the vortices become "heavy", at K = K c they loose sta- 
bility because they attract each other. They exchange a fluxon 
(but nothing is emitted!) and the state turns into a comple- 
mentary one, i.e., \\— *\[. This can be seen in Fig.[2jd): for 



K = 1 .371 the system is still in the J 1 state, but at K = 1 .4jt the 
system is already in the \ [ state. Further increase of K results 
in a decrease of the vortex amplitudes and fluxes as it should 
be for the \ [ state at K — > 27C. The value of K ( depends on a 
and the K c (a) dependence is presented and discussed below. 

B. Where is the flat phase state? 

In the case of 7C-vortices the ground state depends on 
a^^l^m. For a < a c = k/2 the ground state is the so-called 
flat phase state /j = without magnetic flux. For a > a c the 
ground state is made of two AFM ordered semifluxons fj,. In 
addition there is a FM state ||, which exists and is stable for 
any a, but its energy is larger than the one of the AFM or flat 
phase state. 

In the case of arbitrary K, there is no flat phase state at all, 
because /j — 0, in general, is not a solution of Eq. 0, with Q(x) 
given by Eq. d!5t or Eq. d!6t . The only exception is K = Tin, 
with integer n. 

It is therefore interesting to see what happens with the two 
AFM states discussed above when the distance a between 
them becomes smaller than a c . We note here that the value 
of a c =k/2 is valid only for K = 7C. For other K the value of 
a c may change, which is indeed the case for asymmetric AFM 
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n — i — i — i — i — i — i — i — i — i — i — i — t — i — i — i — i — i — i 

1 " 4.CL J— 




discontinuity kIk 

FIG. 3: The magnetic flux A(|) in the right half of the LJJ (x > 0) in 
the symmetric AFM state as a function of K for different values of 
a shown next to each curve. Dashed line A<|) = K shows the limit of 
very large a. 



state, as we will see below. 



1. Symmetric AFM state. 

The crossover distance a\ = a c =7t/2is the same and does 
not depend on K. Therefore below we just use a c to denote it. 
The phase or magnetic field profiles for the symmetric AFM 
state look qualitatively the same for a < a c as well as for a > 
a c for any given K ^ 71, see Fig.|2jd). The only difference is 
the dependence of the amplitude of magnetic field (or of the 
flux) of each vortex vs. K. In Fig.|3]we show the dependence 
of magnetic flux 



A(j)= / fj.x( x )dx=/j(°°)-!j(0) 
Jo 



(18) 



in one half of the junction (in the whole LJJ the flux integrates 
to zero for a symmetric AFM state) as a function of K for 
different distances a. In normalized units adopted here one 
flux quantum is equal to 2k. To plot each curve in Fig. [3] we 
were sweeping K from to 2n and back to 0. The dependence 
of the vortex's maximum magnetic field on K is qualitatively 
similar to the behavior of the flux A(|)(k). 

As K grows starting from 0, the flux associated which one 
vortex grows more or less proportional to K. If a were very 
large (not interacting vortices) then the flux of one vortex 
would be just A(j) = K, as shown by dashed lines in Fig. [3] For 
finite a, the tails of vortices overlap and the fluxes partially 
cancel each other resulting in a smaller slope of A(|)(k). 

For a < a c , the flux reaches its maximum at some value 
< K„, < 71 and starts decreasing (although K increases!), 
smoothly reaching zero at K = 7C. At this point the state of 



the system is the flat phase state /j = 0, which is allowed at 
K = 7t. Further increase of K makes the flux negative. This 
corresponds to the smooth continuous transition from the J ] 
to the I" [ state as K passes the point K = 7C. At K — * 2n the flux 
vanishes as it should be for the \ [ state. When K is swept back 
from 2k to the flux goes back along the same curve without 
hysteresis, again making the \ flat phase state — >J } transi- 
tion at K = 7t. 

For a > a c the flux increases almost all the way up to , 
passing a maximum just before Ke > 7C. At kJ the state 1 J 
becomes unstable and abruptly switches to the complemen- 
tary state \ [ as we already saw in Fig. |2jd). At K — * 2n the 
flux vanishes as it should be for the \ [ state. When K is swept 

M 11 

back from 2% to the state \ [ persists down to K c = 2% — K c , 
as it should be for the complementary state, and the system 
switches back to the original state ] J . 
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FIG. 4: The phase fy(x), magnetic field fi x (x) an d the supercurrent of 
the collective state at a = 1 for different values of K. 

Thus, the crossover which we discovered here is a gen- 
eralization of the crossover between the AFM state and 
the flat phase state which was previously investigated for 
semifluxonsi2i2ii2i The crossover distance a c = n/2 is still 
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the same. It is a bifurcation point, such that for a < a c there 
exists only one symmetric AFM state (J ] or [ \) for given K, 
while for a> a c there are two stable symmetric AFM states (J \ 
and [ \) which exist and are stable for k| ^ < K < TQ For a > a c 
and K outside this interval, again there is only one stable sym- 
metric AFM state: J 1 for < K < k| 1 or f [ for KJ < K < 271. 



2. Asymmetric AFM state 



C. What happens at k^ ? 

The values of Ky (a) and kP (a) are calculated numerically 
and are shown as symbols in Fig. [5] The boundary of ex- 
istence for each solution can be derived using a phase plane 
analysis^. In fact, instead of searching for the critical value 

of K at given a, in the phase plane analysis we search for the 

11 11 

switching distances a s and a s ' as a function of the disconti- 
nuity K. 

First, we treat symmetric FM state. At the edge of the solu- 
tion existence the phase trajectories 



The phase or magnetic field profiles for the asymmetric 

1 I 

AFM state look qualitatively different for a < a c and for 



a > a c . Note, that in this case a c 1 is a weak function of K such 



11 



that a c changes from Jt/2 : 



1 .57 at k = 7C towards a\ ~ 1 



at K — > and K — ► 2k. In Fig.|4]we show the shape of two vor- 
tices in this state for a = 1 < ay and different K [c.f., Fig.|2ja) 
with a > a c ]. We see that at K = 7C the state indeed degener- 
ates into a flux-less flat phase state. For K ^ K, magnetic flux 
appears, but it is distributed symmetrically between the vor- 
tices so that Hx{— x) = /Jx(x). In fact, two vortices behave as a 
single object with the maximum or minimum of the magnetic 
field at x = 0. We call this state a collective state. When K — > 
or K — > 2k the shape of the vortex approaches the shape of a 
single fluxon (antifluxon), again centered at x = 0. In contrast, 

1 L 

for a> a c the asymmetric AFM state at K — > or K — ► 2k re- 
sults in a fluxon (antifluxon) centered at x = ±a/2. We do 
not show the dependence of A0(k) similar to the one shown 
in Fig. [3] (for the symmetric AFM state) because for the 1 [ 
state it is trivial: A0(k) = 2k — 2k, i.e., it is just a straight line 
crossing zero at K = Jt. The flux above refers to the flux in the 
whole LJJ. We note that this dependence follows from simple 
topological considerations and is valid for any a. 

The crossover distance al^K) separates the weakly cou- 
pled vortices (a > a\ ) which are in the asymmetric AFM state 

1 1 

from the collective strongly coupled vortex state (a < a c l ) in 
which the flux is distributed equally and symmetrically be- 
tween the discontinuities. The dependence ai^K) can be cal- 
culated as an edge of existence of the symmetric solution (so- 
lution with maximum a = a^ x for given k). The phase plane 
analysis 23 shows that 



fl H(*) 



(k)=2F -|l-sin- 



(19) 



A) 



= ±y/l -cos(p); 



= i^/C-cos^-K); 
/ii 2K) 0u) = ± v /1-cos( jU -2k), 



(20) 
(21) 
(22) 



corresponding to the 0, K, and 2k regions, touch each 
other. This happens at value of integration constant C = 
1 + 2sin(K/2), and the phase /j in the middle K-region changes 
from (k + tc)/2 at x = -a /2 to K at x — and to (3k — k)/2. 
At this state the value of a can be calculated as 



ay(K) = 2y/m FM 



,Jt. \ / 3tc — K. 
F ( - \m FM J - F I — - — \m FM 



(23) 



with modulus 



m FM (K) 



1 



sin: 



Note, that F(7t/2|m) = K{m) which is a complete elliptic inte- 
gral of the first kind—. The edge of solution existence given by 
Eq. (I23> is plotted in Fig.|5ja) as continuous line. One can see 
the perfect agreement between Eq. d23i and the result of direct 



numerical simulation of Eq. (|2). At a — > 0, k| (a) w 7C + 2a. 
This limiting behavior is clear, as two semifluxons with a = 
form a fluxon. As soon as we try to increase K the state be- 
comes unstable, because a vortex with the topological charge 
larger than 27C is unstable^. 

Second, we analyze the symmetric AFM state. If one draws 
the trajectories 



(2K), s 



±y/i -cos(p); 
±y/C — cos(fi-K); 



= ±0 -cos(^), 



(24) 
(25) 
(26) 



is reached when the phase trajectory corresponding to the mid- 
dle region \x < a/2\ of LJJ is tangential to the phase trajecto- 
ries of the left and the right regions. This happens when the 
phase in the middle region changes from fj(— a/2) = (k — k)/2 
to fi(+a/2) = 3(7C-K)/2. The function F(cp|m) in Eq. dT9l 
is the incomplete elliptic integral of the first kind 24 . Note, the 



limiting behavior a c l {K) = k/2, while a c l (0) = ay (2%) = 
21n(l + V2) w 1.7628, which is in good agreement with nu- 
merical results. 



corresponding to the O-K-0 LJJ on the phase plane, one can see 
that the symmetric AFM solution exists when the integration 
constant C changes from C = 1 (corresponding to a = °o) down 
to C c = 1 — 2sin(K/2) (some finite a) corresponding to phase 
trajectories for and K regions touching each other. The naive 
conclusion that the state with C = C c will correspond to the 
switching (minimum) distance a s is wrong! As it was pointed 
out by Susanto 2 ^, the dependence a(C) is not monotonous and 

has a minimum on the interval C = C c (k) . . . 1. It is this mini- 

i 

mum distance, which corresponds to a s ■ 
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0123456 0123456 
distance between vortices, a distance between vortices, a 



FIG. 5: The behavior of Kc(fl) for symmetric FM (a) and AFM (b) states of two vortices calculated using direct numerical simulation of Eq. 
shown by symbols. The continuous lines show the boundaries of existence regions a] (k) given by Eq. 1231 (a) and a]' (k) (b) calculated using 
the phase plane analysis. The area of existence and stability of each state is shadowed. 



For given K and C the distance a is given by 
fl(C,K) =2 



K| i±£VF(W)|i±^ 



(27) 



when 



.25 



9o(C,k) = sgn(9o(C,K))arcsin 



1 -cos8 (C,k) 



1+C 



and 



K 1 — C 

O (C, K) = 7t - - - arcsin — — ^ 
2 2 sin j 

The value of C at which a(C,K.) has a minimum for given 
K can be found from da(C,K)/dC = 0. One ends up with a 
rather bulky transcendental equation for C. For the sake of 
simplicity we find the minimum value a m ; n (K) = a s (k) for 
given K numerically. The result is shown in Fig. 0b) by a 
continuous line. One can see the perfect agreement between 
the result obtained from the phase plane analysis and the result 
of direct numerical simulation of Eq. (|2j. Note, that kJ — ► Jt 
when a — > a c = Jt/2. 

At large a — > °° both symmetric AFM and FM states be- 
come unstable when each vortex grows up to a fluxon size 
(k = 2;t). This is a natural limit for non-interacting vortices. 



D. Lowest energy state 

One may ask, what is the lowest energy state among the 
four irreducible states discussed above? As we know, for a 
semifiuxon the AFM state (for a > a c ) or the fiat phase state 
(for a < a c ) always has a lower energy than the FM state 19 , 
i.e., for any a. Is this the case for arbitrary vortices, too? 
Before answering this question, we note that the competition 
of energies can only be considered among the states with the 
same structure of discontinuities. 



Thus the first question is: given two discontinuities 
(— K, — k) at a distance a between them, which state has lower 
energy: the asymmetric AFM or the symmetric FM? Look- 
ing at Fig. |2ja)-(b), one can qualitatively say, that at small 
K the symmetric FM state has lower energy than the asym- 
metric AFM state. Instead, at large K the AFM state should 
have lower energy. The numerically calculated energy E(k) 
of both states is plotted in Fig. |6ja) for different values of a. 
To calculate the energy we have solved numerically Eq. (|2j to 
get the solution /j(x) and calculated the energy by integrating 
Eq. (|6) over the junction length. As one can see in Fig.|6la), 
the asymmetric AFM state has lower energy than both sym- 
metric FM states only in a rather narrow interval of K around 
K = 71. If K is outside of this interval one of the symmetric 
FM states has lower and another higher energy than the cor- 
responding asymmetric AFM state at the same K. In Fig. |6ja) 
one can also see transitions at K c corresponding to the switch- 
ing from the "heavy" symmetric FM state to the asymmetric 
AFM state. 

The second question is: given two discontinuities (+K, — k) 
and the distance a between them, which state has lower en- 
ergy, the symmetric AFM or the asymmetric FM? Looking 
at Fig.|2jc)-(d), one can qualitatively say, that at small K the 
symmetric AFM state has lower energy. At larger K it is not so 
clear. The calculations give the E(k) plots shown in Fig.|6jb). 
As we see, in the absence of hysteresis (a < a c ) the symmetric 
AFM state always has lower energy than the asymmetric FM 
state. For large a, the most "heavy" of the symmetric AFM 
states may have larger energy than the corresponding asym- 
metric FM state at the same K. 



V. CONCLUSIONS 

We investigated possible ground states of fractional vortices 
formed at one and two K-discontinuities of the phase. 

In case of one — K-discontinuity we derived the shape of a 




FIG. 6: Numerically calculated energy E(k) of competing states for different values of a, shown next to each curve, (a) asymmetric AFM vs. 
symmetric FM; (b) symmetric AFM vs. asymmetric FM; 



direct +K-vortex, see Eqs. 10 and dlOi . and a complementary 
(k — 27c)-vortex, see Eqs. (II It and d!2i . as well as their energy 
as a function of K dl4> . In the general case K ^ K these vortices 
are not mirror symmetric like semifluxons. 

Due to such a broken symmetry, the ground state of the sys- 
tem with two discontinuities consists of four different states: 
asymmetric and symmetric AFM and FM states. Symmet- 
ric states consist of two direct vortices or two complemen- 
tary vortices, while asymmetric states consist of one direct 
and one complementary vortex. Asymmetric states are sta- 
ble for < K < 271 and their properties are symmetric (en- 
ergy is the same, magnetic field just changes sign) with re- 
spect to the transformation K — > 2k — K. Symmetric states 
exist only in some interval of K. For example, the state J] 

I 1 

exists only for < K < K c > n and the state 1 ] exists only 
for < k < kJ 1 > 7C. The values (a) and kJ 1 (a) have been 
found, see Fig. |5] If K exceeds K c , the state turns itself into 
another state, e.g., ] J— > [ \ or ]] — >1 [. For the symmetric AFM 
state a new, more general, meaning of the crossover distance 
Oc = a c = tc/2 between the discontinuity points is discov- 
ered. If a < a c , the transition between the ] J and [ \ states is 
smooth, otherwise it is abrupt and is associated with the insta- 
bility of the ] \ state at K = 1Q . The transitions between differ- 
ent weakly coupled states which can be induced by changing 
K are summarized in Fig.Q 

The asymmetric AFM state exists only for a > a c (k), 
while for a < a\ (k) it turns itself into a new strongly cou- 
pled "collective" state which is perfectly symmetric and can 
not be simply constructed out of the single vortex states. 
The crossover distance a c is a weak function of K given by 
Eq. (El. 

Finally, we have calculated and compared the energies of 



competing states. We showed that the energy of the asymmet- 
ric AFM state can be larger as well as smaller than the one of 
the symmetric FM state depending on K. For fixed K it may 
or may not depend on a. On competition between symmetric 
AFM and asymmetric FM states, there is always (for any a 
and k) a symmetric AFM state which has lower energy than 
the asymmetric FM state. The details are presented in Fig. [6] 

As we see, the variety of the ground states of arbitrary vor- 
tices is much more rich than the one of semifluxons. The 
asymmetry may be exploited in information processing de- 
vices based on distinction and controllable manipulation of 
the vortex states. 

The case of two vortices represents the simplest system of 
two coupled vortices. Of course, in the case of N > 2 discon- 
tinuities, one should expect even a larger variety of possible 
vortex states which may strongly depend on the parity of N, 
but the treatment of those configurations should depend on 
practical needs. In this sense the case of two vortices is espe- 
cially important as it is the first candidate for implementation 
of classical or quantum bits based on fractional vortices. The 
difficulty with one vortex is that it emits a fluxon when it flips. 
The AFM state, instead can flip like 1 J ] or J ["«-► \\ without 
any emission. 
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FIG. 7: Possible transitions between different weakly coupled two vortex states which can be achieved by varying K. 
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